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Abstract 

The indication for the a decay of 180 W with a half-life T 1 Q /2 =l.ll|]j(stat)±0.3(syst)xl0 18 yr has 
been observed for the first time with the help of the super-low background 116 CdWO"4 crystal 
scintillators. In conservative approach the lower limit on half-life of 180 W has been established 
as T 1 " 2 ( 180 W) > 0.7xl0 18 yr at 90% C.L. Besides, new Ty 2 bounds were set for a decay of 

182 W) 183 W) 184 W and 186 W &t the levd Q f 1Q 20 yr _ 

PACS number(s): 23.60.+e; 29.40.Mc 
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1 INTRODUCTION 



Alpha decay is allowed energetically for the five naturally occurring isotopes of tungsten [1], 
but it was never observed up to now. The a activity of tungsten with a particle energy of 
about 3 MeV and half-life T\/2 = 2.2x5xl0 17 yr (where 5 is the relative abundance of isotope) 
was declared in an early experiment with nuclear emulsion technique [2]. Because Q a value 
for all natural tungsten nuclides, except 180 W, is less than 2 MeV, the result of ref. [2] could 
be attributed to a decay of 180 W (5 « 0.12% [3]) with T 1/2 = 2.6xl0 14 yr. However, this 
observation was ruled out in the work [4], where the cadmium tungstate (CdW0 4 ) crystal 
scintillator (mass of 20.9 g) was used as a source and detector of a decay events simultaneously. 
After 193 h of measurements the limit T\ji > l.OxlO 15 yr was established [4]. A similar 
restriction (T 1/2 > 9.0xl0 14 yr) was also obtained in experiment (66.7 h of exposition) with 
ionization counter (1200 cm 2 in area) and thin (83 /ig/cm 2 ) sample of W (with total mass of 
79 mg) enriched in 180 W to 6.95% [5]. 

These bounds were improved only recently in the measurements with two scintillators: 
CdWO-4 (mass of 452 g, running time of 433 h), and 116 CdW0 4 enriched in 116 Cd to 83% (91.5 
g, 951 h) [6]. The limits on the half-lives for a decay of different W isotopes were set in the 
range of ~ 10 17 -10 19 yr (see Table). 

In this paper the new results of Kiev-Firenze experiment (2975 h of exposition) on a decay of 
natural tungsten isotopes are described (the preliminary analysis of the 975 h data was presented 
in [9]). They were obtained with the help of the super- low background spectrometer [10] based 
on enriched 116 CdW04 crystal scintillators. The sensitivity of this apparatus to measure a 
activity was enhanced substantially in comparison with our previous work [6] mainly due to 
the developed technique of pulse-shape analysis of the data, which allows us to distinguish 
events caused by a particles from those by 7 rays (/3 particles). 

1 Current address: INFN, Laboratori Nazionali del Gran Sasso, 67010 Assergi (AQ), Italy 
2 Corresponding author. E-mail address: zdesenko@kinr.kiev.ua 
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Table 1: Theoretical and experimental half-lives (or limits at 90% CL.) for a decay of natural 
W isotopes obtained in present work. The uncertainties of calculated T1/2 values are related 
with uncertainties of the Q a . The most stringent previous experimental bounds [6] are cited 
for comparison. 



Isotope, 
abundance 
[3] 


Q a , MeV 
[1] 


Calculated 7i/2, yr 


Experimental 7\/2, yr 


based on [7] 


based on [8] 


present work 


[6] 


180 W 
0.12(1)% 


2.516(5) 


8.3ii;«xl0 17 


2.0i^xl0 18 


l.ii^xio 18 

> 0.7xl0 18 


> 7.4xl0 16 


182 W 

26.50(16)% 


1.774(3) 


3.0 xlO' 2 


l.^xlO 33 


> 1.7xl0 20 


> 8.3xl0 18 


183 W 

14.31(4)% 


1.682(3) 




> 5.7t 1 ; 3 ,xl0 38 


> 0.8xl0 20 


> 1.9xl0 18 


30.64(2)% 


1.659(3) 


3.8^ xlO 35 


2.6t^xl0 36 


> 1.8xl0 20 


> 4.0xl0 18 


186 W 

28.43(19)% 


1.123(7) 


8.7^ xlO 55 


2.o^;5 xio 57 


> 1.7xl0 20 


> 6.5xl0 18 



2 EXPERIMENT AND DATA ANALYSIS 

2.1 Experimental set-up 

The set-up, used in our study, was originally devoted to the search for the neutrinoless double 
beta decay of 116 Cd [10]. The 2(3 decay experiment is carried out by the INR (Kiev) 3 in the 
Solotvina Underground Laboratory (allocated in a salt mine 430 m underground [11]) and 
results were published elsewhere [10, 12]. The high purity 116 CdW04 crystal scintillators, 
enriched in 116 Cd to 83%, were developed and grown for the search [12]. Their light output 
is ~30-35% of Nal(Tl). The fluorescence peak emission is at 480 nm with principal decay 
time of ?sl4 /is [13]. The CdW0 4 refractive index equals 2.3. The density of crystal is 7.9 
g/cm 3 , and material is non-hygroscopic and chemically inert. In the apparatus (see for details 
[10]) four 116 CdW04 crystals (total mass 330 g) are exploited. They are viewed by a low 
background 5 EMI phototube (PMT) with RbCs photocathode through light-guide 010x55 
cm, which is glued of two parts: quartz (25 cm) and plastic scintillator (Bicron BC-412, 30 
cm). The enriched 116 CdW0 4 crystals are surrounded by an active shield made of 15 natural 
CdW0 4 crystals of large volume with total mass of 20.6 kg. These are viewed by a PMT 
through an active plastic light-guide 017x49 cm. The whole CdW04 array is situated within 
an additional active shield made of plastic scintillator 40x40x95 cm, thus, together with both 
active light-guides, a complete Anx active shield of the main ( 116 CdW0 4 ) detector is provided. 

The outer passive shield consists of high purity copper (3-6 cm), lead (22.5 - 30 cm) and 
polyethylene (16 cm). Two plastic scintillators (120x130x3 cm) are installed above the passive 
shield and are used as cosmic muon veto counters. The set-up is carefully isolated against 
environment radon penetration. All materials for the installation were previously tested and 

3 From 1998 by the Kiev-Firenze collaboration [10]. 



2 



selected for low radioactive impurities in order to reduce their contributions to background. 

An event-by-event data acquisition system records the amplitude, arrival time, additional 
tags (the coincidence between different detectors) and pulse shape (in 2048 channels with 50 
ns channel's width) of the 116 CdW0 4 detector in the energy range of 0.08 — 5 MeV. 

The energy scale and resolution of the spectrometer were determined with the 7 sources 
22 Na, 40 K, 60 Co, 137 Cs, 207 Bi, 226 Ra, 232 Th and 241 Am. The energy dependence of the resolution 
in the energy interval 60-2615 keV is expressed as follows: FWHM 7 = -44+^2800 + 23.4 • £ 7 , 
where energy E 1 and FWHM 7 are in keV. For instance, energy resolutions of 33.7%, 13.5%, 
11.5% and 8.0% were measured for 7 lines with the energies of 60, 662, 1064 and 2615 keV, 
respectively. The routine calibration was carried out with a 207 Bi (weekly) and 232 Th (monthly) 
7 sources. The dead time of the spectrometer and data acquisition was permanently controlled 
with the help of a light emitting diode optically connected to the main PMT (typical value was 
about 14%). 

Due to active and passive shields, and as a result of the time-amplitude [10] and pulse-shape 
analysis of the data [13], the background rate of 116 CdW0 4 detectors in the energy region 2.5 
- 3.2 MeV (0i/2/3 decay energy of 116 Cd is 2.8 MeV) was reduced to 0.04 counts/yr-kg-keV. It 
is the lowest background which has ever been reached with crystal scintillators. 

2.2 Response of the 116 CdW04 detector to a particles 

It is well known that relative ratio of the CdW0 4 scintillation light output for a and (3 particles 
with the same energies (so-called a/f3 ratio 4 ) depends on the energy of the absorbed particles 
[6]. Because we are looking for the a decay of W isotopes, such a dependence must be precisely 
measured in the a energy range of 1-3 MeV (see Table). Unfortunately, among nuclides from 
U/Th families there are no a emitters with such energies (the lowest available a energy is 4.0 
MeV from 232 Th). To overcome this problem, a special method of calibration was developed, in 
which a collimated beam of a particles from 241 Am source was passed through a thin absorber 
with known thickness, thus the energy of a particles after absorber can be calculated precisely. 
Furthermore, it was measured with the help of a surface-barrier semiconductor detector. The 
dimensions of the collimator are 00.75 x 2 mm, and thickness of a single mylar film absorber is 
0.65 mg/cm 2 . By using this technique and different sets of absorbers, a particles with energies 
of 0.46, 2.07, 3.04, 3.88, 4.58, and 5.25 MeV were obtained, which allows us to calibrate our 
detector in the energy range of interest. 

It is also known that the light output of crystal scintillators may depend on the direction 
of a irradiation relative to the crystal axes [14]. To study this effect for CdW0 4 , 116 CdW0 4 
and CaW0 4 crystals, they were irradiated by a particles in three directions perpendicular to 
the (010), (001) and (100) crystal planes 5 , and three experimental dependences of the a/f3 
ratio (corresponding to three mentioned directions) as a function of a particle energy were 
derived from measurements. However, in a real crystal the amplitude of a light signal depends 
also on the point, from which scintillation light is emitted (so-called nonuniformity of light 

4 The detector energy scale is measured with 7 sources, thus the notation "a/7 ratio" could be more adequate. 
However, because 7 rays interact with matter by means of the energy transfer to electrons, in present paper we 
are using traditional notation "a//3 ratio". 

5 In the 116 CdWC>4 crystal (010) plane is perpendicular to the cylinder axis, so for the (010) direction a source 
was placed on the top of crystal (in the center of the flat circle surface). For the (001) and (100) directions 
crystal was irradiated in the middle of cylindrical surface. 
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collection). The latter can distort the anticipated effect of crystal's orientation, and hence, 
should be properly taken into account. With this aim, the light propagation in the CaW0 4 
scintillator, for which the effect of crystal orientation was not observed, was simulated with 
the help of the GEANT3.21 package [15] for the light emitting points on the top and side 
surfaces of the crystals (with dimensions 20x20x11 and 40x34x23 mm), as well as for those 
uniformly distributed inside the crystal. It was found that results of simulations are in a good 
agreement with experiment. Then, the same simulations were performed for the 116 CdW04 
crystal (032 x 19 mm), and results of each calibration measurement with this crystal were 
corrected by using simulated distributions of the light collection for a particles emitted from 
the corresponding point. So, values of the a//3 ratio for the direction 1 was multiplied by 
a factor of 0.85, while for the direction 2 and 3 - by a factor of 1.08. As an example, such 
corrected dependences of the a/ f3 ratio versus the energy and direction of a particles are shown 
in Fig. 1 for enriched 116 CdW0 4 crystal (032 x 19 mm). 




I I I | I I I | I I I | I I I | I I I 

2000 4000 6000 8000 10000 

Energy of a particles (keV) 

Figure 1: The dependence of the a/ (3 ratio on energy and direction of incident a particles 
measured with the enriched 116 CdW04 scintillator (032 x 19 mm). The crystal was irradiated 
by external a source in directions perpendicular to (010), (001) and (100) crystal planes (dir. 
1, 2 and 3, respectively). Besides, internal a peaks of 224 Ra, 220 Rn, 216 Po and 212 Po from the 
intrinsic contamination of the 116 CdW0 4 crystal were used (see text). Solid lines represent the 
fit of ct/P ratio dependence. In the insert the behaviour of a/ (3 ratio measured with CdW0 4 
scintillator (025 x 0.9 mm) in the direction 1 is shown. 

In addition to the measurements with the external source, the a peaks from the intrinsic 
trace contamination of the 116 CdW0 4 crystals by nuclides from Th chain were also used for 
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determination of the a//3 ratio. These peaks were selected from the background by using the 
time-amplitude analysis [10, 16]. For instance, the following sequence of a decays from the 
232 Th family was searched for and observed: 224 Ra (Q a = 5.79 MeV, T 1/2 = 3.66 d) -> 220 Rn 
(Q a = 6.40 MeV, T 1/2 = 55.6 s) -> 216 Po (Q a = 6.91 MeV, T 1/2 = 0.145 s) -> 212 Pb. The 
obtained a peaks (a nature of events was confirmed by a pulse-shape analysis described below), 
as well as the distributions of the time intervals between events are in a good agreement with 
those expected for a particles of 224 Ra, 220 Rn and 216 Po (see Fig. 2). On this basis the activity 
of 228 Th in 116 CdW0 4 crystals was determined as 39(2) yuBq/kg 6 . Moreover, the a peak of 
212 Po (the daughter of the 212 Bi) was reconstructed with the help of the front edge analysis of 
the scintillation signals. The energy and time distributions of the sequence of (3 decay of 212 Bi 
and a decay of 212 Po, selected from the background, are presented in Fig. 3 (a, b and c), while 
typical example of such an event is shown in Fig. 3d. 
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Figure 2: The a peaks of 224 Ra, 220 Rn, 216 Po selected by the time-amplitude analysis from 
background data accumulated during 14745 h with the 116 CdW0 4 detector. In the inserts the 
time distributions between the first and second (and between second and third) events together 
with exponential fits are presented. Obtained half-lives of 220 Rn and 216 Po (6lig° s and 0.144(8) 
s, respectively) are in a good agreement with the table values: 55.6(1) s and 0.145(2) s [21]. 



The values of the a//3 ratio, derived with the internal a peaks, and their fit are depicted in 
Fig. 1. Fit yields for the energy range 2.0 - 8.8 MeV: a//3 = 0.083(9) + 0.0168(13) • E a , where 

6 The same technique was applied to the sequence of decays from the 235 U anc l 238y families. Activity of 
5.5(14) ^Bq/kg for 227 Ac ( 235 U family) and limit < 5 [iBq/kg for 226 Ra ( 238 U family) in the 116 CdW0 4 crystals 
were set. 
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Figure 3: The energy (a, b) and time (c) distributions for fast sequence of (3 ( 212 Bi, Q@ = 2254 
keV) and a ( 212 Po, E a = 8784 keV, T1/2 = 0.299(2) /j,s) decays selected from the background 
data by the pulse-shape analysis, (d) Example of such an event in the 116 CdW04 scintillator. 



E a is in MeV. 7 Because in measurements with internal sources the effect of crystal's orientation 
is averaged, such an extrapolation of the fit into the energy region 2.0 - 5.5 MeV is reasonable 
and it was also proved by the behaviour of the a/f3 ratios measured with external sources (see 
Fig. 1). At the same time, in the a energy interval of 0.5 — 2.0 MeV we find that a//3 ratio is 
decreased with the energy. It was also confirmed by the measurement with the thin (025 x 0.9 
mm) CdW0 4 crystal scintillator (see insert in Fig. 1). Similar behaviour (fall with energy in 
the 10 — 100 keV energy region) of relative scintillation efficiency for Ca and F recoil nuclei in 
CaF 2 (Eu) scintillator was reported in [18]. Thus, for 0.5 — 2.0 MeV energy range we obtain: 
a/(3 = 0.15(3) - 0.015(8) • E a , where E a is in MeV. 

The calibration data were also used to determine the energy resolution of the detector for 
a particles: FWHM a (keV)= 33 + 0.247 E^, where is the energy of a particles in 7 scale 
expressed in keV. 

7 The growth of the a/ (3 ratio with a particles energy has been earlier observed for cadmium tungstate [6] 
and other scintillators [14, 16, 17]. 
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2.3 Pulse-shape analysis 

The pulse-shape analysis of CdW04 scintillation signals was developed on the basis of the 
optimal digital filter [19], and clear discrimination between 7 rays (electrons 8 ) and a particles 
was achieved [13]. The pulse shape of cadmium tungstate scintillators can be described as: 
f(t) = S Mjiji — r o) • (e~*/ Ti — e~*/ r °), where Ai are amplitudes and Tj are decay constants for 
different light emission components, r is the integration constant of electronics (?» 0.18 fis). To 
provide an analytic description of the a or 7 signals, the pulse shape resulting from the average 
of a large number of individual events has been fitted with the sum of three (for a particles) or 
two (for 7-s ) exponents, giving the reference pulse shapes f a (t) and / 7 (t) (see for details ref. 
[13]). For the enriched crystals (used in the experiment) the following values were obtained: 
A?=80.9, rf=12.7 /is, A£=13.4, r 2 a =3.3 fis, A%=5.7, r 3 Q =0.96 /is for ^5 MeV a particles and 
A]=94.2, 7-7=13.6 fjs, Al=5.6, r 2 7 =2.1 ^s for «1 MeV 7 quanta. 

In the data processing the digital filter was applied to each experimental signal f(t) with 
aim to obtain the numerical characteristic of its shape defined as SI (shape indicator): SI = 
J2f(tk) x P{tk)/ J2 f(tk), where the sum is over time channels k, starting from the origin of 
signal and up to 50 fis, f(t k ) is the digitized amplitude (at the time t k ) of a given signal. The 
weight function P(t) is defined as: 

^(*) = {7a(*)-7 7 (*)}/{7«(*) + 7 7 (*)}- 

The SI distributions measured with different a and 7 sources (see some examples in 
Fig. 4) are well described by Gaussian functions, whose standard deviations cr(SI a ) and 
a(SI^) depend on energy: cr(SI a ) =-0.3+0.7xl0~ 3 i?2+1700/i?2 for a particles, and a(S7 7 ) = 
1.51— 0.23x 10 _3 -E 7 +402/i? 7 for 7 quanta, where E% and E 1 are in keV. As it is seen from Fig. 
4, electrons (7 rays) and a particles are clearly distinguished for the energies above 0.6 MeV 
(E a ~3 MeV), while pulse shape discrimination ability of CdW0 4 scintillator is worsening at 
lower energies. Nevertheless, it is also visible from Fig. 4a that even 2 MeV a particles (E^ 
paO.3 MeV) can be separated from 7 background with reasonable accuracy. For instance, with 
requirement that ^87% of a particles must be registered we will get about 13% of 7 events as 
background (see Fig. 4a). 

The dependence of the pulse shape on energy and direction of a particles measured with 
the enriched 116 CdW04 crystal scintillator (032x19 mm) is presented in Fig. 5. In the energy 
range of 0.5 — 7.0 MeV the average dependence of the shape indicator on the energy can be 
approximated by the function: SI a = 29.5 — 0.195x 10~ 2 E^, where E^ is in keV. 

For 7 quanta the energy dependence of the shape indicator was measured with 7 sources in 
the energy range of 0.04 - 3.2 MeV: 57 7 = 18.4 - 0.117x 10~ 2 £ 7 + 0.54x lO" 6 ^, where E 1 is 
in keV. 

Besides, a digital filter for the pulses of the plastic scintillator light guide was developed 
and clear separation was obtained for the pure events in the plastic and cadmium tungstate 
scintillators. It allows us to discriminate plastic pulses in more complicated cases, when they 
are overlapped with the signals of 116 CdW0 4 crystals. The use of this filter leads to some loss 
of the registration efficiency for the events in the 116 CdW04 detector (of the order of ~5%), 
which, however, can be correctly taken into account on the basis of calibration measurements. 

8 Bccause 7 rays interact with matter by means of the energy transfer to electron, it was assumed and 
experimentally proved with j3 particles from the decay of internal 113 Cd that pulse shapes for electrons and 7's 
are the same. 
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Figure 4: The examples of the shape indicator distributions measured by the enriched 
116 CdW04 crystal scintillator (032x19 mm) with a particles and 7 rays with the different 
energies: (a) £ 7 = 0.2-0.4 MeV, E a = 2.07 MeV; (b) £ 7 = 0.5-0.7 MeV, E a = 3.04 MeV; (c) 
E y = 0.8-1.0 MeV, E a = 3.88 MeV; (d) £ 7 = 1.2-1.4 MeV, E a = 5.25 MeV. 



3 RESULTS AND DISCUSSION 

3.1 Background interpretation 

The background spectra of (7 + f3) and a events measured by four 116 CdW04 crystals (330 g, 
exposition 2975 h) are depicted in Fig. 6. 

The 7(/9) spectrum shown in Fig. 6 was built by selecting the following interval of SI 
values: S7 7 — 2.4<x 7 < SI < SI-, + 2.4<7 7 , which contains 98% of ^((3) events. In the low energy 
region the background of the 116 CdW0 4 detector is caused mainly by the fourth-forbidden (3 
decay of 113 Cd (T 1/2 = 7.7 x 10 15 yr [20], = 316 keV [21]) which is present in the enriched 
crystal with abundance of 5 ~ 2%. The distribution above ~ 350 keV is described by a trace 
contamination of the enriched and shield crystals by 40 K, 137 Cs, 113m Cd, two neutrino double 
beta decay of 116 Cd with Ti/ 2 = 2.6 x 10 19 yr [10], and external 7 rays. 

The energy spectrum of a particles (Fig. 6) was obtained by selection of events with the 
following values of shape indicator: S7 7 + 4<r 7 < SI < SI a + 2Aa a . Under such restrictions the 
efficiency of the pulse-shape analysis (t)psa) depends on the energy of a particles. For example, 
for a peak of 180 W this efficiency equals 49.5%, while additional use of the filter for the plastic 
pulses discrimination decreases this value down to 47%. 
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Figure 5: The dependence of the shape indicator on the energy and the direction of incident 
a particles measured with the enriched 116 CdW04 crystal scintillator. The a peaks of 224 Ra, 
220 Rn and 216 Po were selected by the time- amplitude analysis of background data. 



Taking into account the fact that secular equilibrium in crystal scintillators is usually broken, 
the distribution of the a events is well reproduced by the model (see insert in Fig. 6), which 
includes a decays from 232 Th and 238 U families. For illustration, the results of the pulse-shape 
analysis of the data (for energy above 350 keV) are presented in Fig. 7 as three-dimensional 
distribution of the background events versus energy and shape indicator. In this plot one can 
see clearly separated population of the a events, which belong to U/Th families. The total 
a activity associated with the peak in the energy region 400 — 1500 keV is 1.40(10) mBq/kg. 
However, because of a poor energy resolution for a particles and uncertainty of the a/ (3 ratio, we 
give only limits on the activity of nuclides from uranium and thorium families in the 116 CdW04 
scintillators received from the fit of the spectrum in the range 400 — 1500 keV: 232 Th < 0.15 
mBq/kg, 238 U ( 234 U) < 0.6 mBq/kg, 230 Th < 0.5 mBq/kg, 210 Pb < 0.4 mBq/kg. 

The low energy part of the a spectrum (below 200 keV) can be explained by the PMT noise, 
residual 7(/3) background, 7 or (3 events in the 116 CdW04 scintillators with small admixture 
of the plastic light guide pulses (which were not discriminated by the pulse-shape analysis), 
decays of Th and U daughter a nuclides located inside defects or inclusions of the crystals (in 
that case alpha particles can lose a part of their energy without scintillation light emission), 
etc. 

Lastly, a small peak is visible in the a spectrum of Fig. 6 at the energy around 300 keV. 
Since the 2.46 MeV a peak of 180 W is expected at 307 ± 24 keV (with FWHM = 110 keV), 
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Figure 6: Energy distributions of 7(/3) and a events, which were selected by the pulse shape 
analysis from the data of 116 CdW04 crystals (330 g) measured during 2975 h. In the insert, a 
spectrum is depicted together with the model, which includes a decays from 232 Th and 238 U 
families. The total a activity in the 116 CdW0 4 crystals is 1.40(10) mBq/kg. 



and because the position of a peak of 232 Th (with the lowest a energy of 4.0 MeV among all a 
emitters from U/Th families) must be at ^600 keV, it is unlikely that this small peak can be 
attributed to origin other than a decay of 180 W. However, because a platinum crucible was used 
for the growth of our crystals, we analyze the possible imitation of the effect by a activity of 
190 Pt (abundance of 190 Pt is 0.014% [3], T 1/2 = 6.5xl0 n yr, E a = 3164(15) keV [21] or 435±30 
keV in 7 scale). Calculation shows that the detected peak can be caused by platinum pollution 
(homogeneously spread in the volume of the 116 CdW04 scintillators) at the level of ~ 3 ppm. 
To estimate the actual platinum contamination in cadmium tungstate crystals, the results of 
previous low background measurements performed by the Milano-Kiev collaboration with the 
CdW0 4 crystal 9 of 58 g [22] were considered. This experiment has been carried out in the Gran 
Sasso Underground Laboratory, and CdWO^ crystal was used as a low temperature (^25 mK) 
bolometer. Energy resolution (FWHM) of the detector was equal to 5 keV at 2.6 MeV. No 
events were registered in the energy region of 3100 — 3300 keV during 340 h, which allows us 
to set bound on the 190 Pt activity. In accordance with [23, 24] it yields 1.3 counts as limit for 
number of events, which can be excluded with 68% C.L., hence one can set a bound for platinum 
concentration in cadmium tungstate crystal as 1.2 ppm. Besides, two samples of CdW0 4 crystal 
with dimensions 1.5x1.5x0.1 cm were searched (with the help of the electron microscope) for 

9 This crystal was produced in the same apparatus as enriched ones. 
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Figure 7: Three-dimensional distribution of the background events (2975 h of exposition with 
the 116 CdW04 crystals) versus energy and shape indicator. The population of a events belong- 
ing to the U/Th families is clearly separated from j(/3) background. 



inclusions, whose elemental composition is different from that of CdW0 4 . With this aim flat 
surface (2 cm 2 ) of each sample was scanned, and if such an inclusion was observed, the electron 
beam was concentrated on it and the emitted X rays were analyzed with a crystal spectrometer 
(energy resolution is better than 0.1 eV) tuned on a characteristic platinum X rays. For any 
observed inclusions (with diameters in the range 2-30 //m) no emission of platinum X rays was 
found. It results in conclusion that limit (95% C.L.) on platinum concentration in CdW0 4 
crystal is lower than 0.1 ppm (or 0.03 ppm) for Pt inclusions with diameters less than 5 /im (or 
3 jum) . Similar analysis performed with the small sample (~3 cm in diameter) of the enriched 
116 CdW0 4 crystal gives even more stringent restrictions 10 : 0.07 ppm (or 0.02 ppm). 

Therefore, one can conclude that all bounds obtained for homogeneous Pt contaminations 
in the CdW0 4 crystals are well below the level of 3 ppm, at which a peak of 180 W could be 
imitated. Moreover, it is also known that in crystals, which were grown in platinum crucible, Pt 
is not distributed homogeneously, but it is present in form of precipitates with size around 20 
(Ltm [25] , hence a broad energy distribution instead of an a peak would be observed. The latter 
was proved by our Monte Carlo simulation of the 190 Pt a decays in CdW0 4 crystal with the 
help of the GEANT and event generator DECAY4 [26]. It gives that effect of 180 W a activity 
could be imitated by 190 Pt alpha decays only in the case of platinum particles of 2 — 3 /im size, 
with Pt average concentration in the crystals at the level of 4 — 6 ppm, which is much larger 

10 Due to larger scanned area and better spatial resolution of the electron microscope used in this case. 
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than our experimental limits. 

Nevertheless, because it is impossible to exclude (at least in principle) some other explana- 
tions of the 300 keV peak in the background spectrum, we can treat our experimental result 
only as the first indication for the possible a decay of 180 W. Obviously, final confirmation of 
its existence would be obtained with the help of CdWC>4 crystals enriched/depleted in 180 W. 
However, it is clear that preparation and performance of such a measurement would require a 
strong additional efforts and perhaps a long time. 

3.2 Alpha activity of 180 W and other tungsten isotopes 

Assuming that observed peak at the energy around 300 keV is really caused by a decay of 
180 W, let us estimate its half-life. With this aim the experimental a spectrum in the energy 
region of interest was fitted by a Gaussian distribution (FWHM =110 keV), which represents 
the effect, and by the background model. The latter was built up as a sum of the a peak of 
232 Th and exponential function. Position and area of the searched peak, the thorium a peak 
area, and constants of the exponent were taken as free parameters of the fit procedure, which 
was performed in the energy region (140±20) — (510±30) keV. The best fit (x 2 /n.d.f. = 0.52) 
achieved in the interval 140 — 520 keV (see Fig. 8) gives the area of the searched peak S = 64±26 
counts and its position at 326±15 keV. Therefore, measured energy of 180 W alpha particles is 
2580 ± 290 keV, which is in a reasonable agreement with the theoretical value -E a =2460 keV. 
Taking into account number of 180 W nuclei in the crystals (6.57xl0 20 ) and the total efficiency 
(47%), we get the corresponding half-life value: Tf /2 ( 180 W) = l.lj£!(stat)±0.3(syst) x 10 18 yr. 
The systematic error is related mainly with an uncertainty of the background model. 

In more conservative approach we use results of our fit in order to set the upper limit on 
half-life of 180 W: lim Tf /2 ( 180 W) > 0.7xl0 18 yr at 90% C.L. 11 

In addition, since there are no structural features in the experimental a spectrum, which 
could indicate an a activity of other tungsten isotopes, half-life bounds for these processes were 
estimated. The numbers of candidate nuclei in the detector with mass of 330 g are: 182 W - 
1.45xl0 23 , 183 W - 7.83xl0 22 , 184 W - 1.68xl0 23 , 186 W - 1.56xl0 23 . The position of a peak in 
the 7 equivalent scale, the expected peak width and the efficiency of the pulse-shape analysis 
are: 182 W (E% = 215±57 keV; FWHM a = 86 keV; rj PSA = 35%), 183 W {El = 206±54 keV; 
FWHM a = 83 keV; 7] PSA = 34%), 184 W (E% =204±53 keV; FWHM a = 83 keV; r] PSA = 34%), 
186 W {El = 147±34 keV; FWHM a = 69 keV; t] PS a = 28%). Fitting the experimental a 
spectrum by sum of expected peak and background model (exponential function plus a peaks 
of 180 W and 232 Th) we get the following half-life limits (at 90% C.L.): 

Ti /2 ( 182 W) > 1.7xl0 20 yr, T 1/2 ( 183 W) > 0.8xl0 20 yr, 

Ti /2 ( 184 W) > 1.8xl0 20 yr, T 1/2 ( 186 W) > 1.7xl0 20 yr. 

The obtained experimental results for a activity of naturally occurring tungsten isotopes 
are summarized in Table. 

11 Similar bound T", 2 > lxlO 18 yr at 90% C.L. can be extracted from the mentioned measurements (340 h) 
with the CdWC-4 crystal of 58 g used as a low temperature bolometer [22] , where no events were observed within 
the energy interval 2516±30 keV. 
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Figure 8: Fragment of the a spectrum measured with 116 CdW04 detector during 2975 h together 
with the fitting curve (solid line). The a peak of 180 W with area of 64 counts corresponds to 
half-life 1.1 xlO 18 yr. 



3.3 Comparison with theory 

To our knowledge, there is only one theoretical work, based on the microscopical approach, in 
which the half-life for a decay of 180 W was calculated [27] . It takes into account the systematic 
behaviour of the reduced a widths of even-even nuclei with numbers of neutrons and protons 
N > 84, Z < 84 and penetration factors obtained from realistic cluster wave functions. The 
derived result is: Ty 2 = 7.5 xlO 17 yr and uncertainty of this value is estimated as less than a 
factor of 3 [27]. For other tungsten isotopes theoretical predictions are absent. 

We have calculated the half-lives of all W isotopes for a decay, first, with the help of the 
cluster model [7], which was very successful in describing the a decays of even-even nuclei. For 
example, it reproduces the Ty 2 values mainly within a factor of 2 for a wide range of nuclides 
(from ^Te to iQgHs) and for T", 2 from 10~ 7 s to 10 16 yr. Because a decays of natural W 
isotopes (except 183 W) occur without changes in nuclear spin and parity, we have chosen a set 
of 17 nuclides (from ^Nd to ^Cm) with Ty 2 > 10 5 yr [28], whose a decays also proceed 
without changes in nuclear spin and parity. While comparing a theory with an experiment, it 
is practical to use the value of the deviation x — max(_R, l/R), where R = Ty 2 /Ty 2 . For the 
chosen set of nuclides the cluster model [7] gives a quite reasonable average value of x — 1-9. 
Then, half-lives for W isotopes were calculated with the help of the cluster model [7] and 
obtained results are listed in Table (given uncertainties are related only with uncertainties of 
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Q a ). In particular, T 1 a /2 ( 180 W) = 8.3xl0 17 yr, that is very close to the value of 7.5xl0 17 yr 
derived in [27]. 

Semiempirical relationships are often more effective in the prediction of T^ 2 than micro- 
scopically based calculations. We found in the literature [8, 29, 30, 31, 32, 33, 34, 35, 36, 37] 
18 semiempirical formulae which can be used for our purpose. All of them were also tested 
with the same set of 17 nuclides with T", 2 > 10 5 yr. The best result (average deviation of the 
calculated values from experiment x = 1-9) was found for the relationship of [8] based on phe- 
nomenological fission theory of a decay and valid not only for even-even but also for odd-even, 
even-odd and odd-odd nuclei. The values of Ty 2 obtained in such approach for W isotopes 
are given in Table, where uncertainties are caused, as in the previous case, by uncertainties of 
Q a . For the 183 W decay it was also taken into account that change in parity will additionally 
suppress the decay rate as compared with that of [8]. For 180 W our result Ty 2 = 2.0xl0 18 yr 
is also close to that of [27]. 

Thus, we can conclude that half-life value (T^ 2 = I.IIq^xIO 18 yr) for possible a decay 
of 180 W (indication for which is observed in present work) is in a good agreement with the 
theoretical predictions: Ty 2 = 0.75xl0 18 yr (microscopic approach [27]), and calculated here 
on the basis of semiempirical formula [8] and cluster model [7] as Ty 2 = 2.0 xlO 18 yr and 
Ty 2 — 0.83xl0 18 yr, correspondingly. 

4 CONCLUSIONS 

In the present work the pulse shape and a/ (3 ratio of the CdW0 4 crystal scintillators have been 
studied for three directions of the collimated beam of a particles relatively to the main crystal 
axes in the energy range of 0.5 — 5.2 MeV. The dependences of a//3 ratio and pulse shape on 
the direction of a irradiation have been found and measured. 

By using the super-low background 116 CdW04 crystal scintillators, the indication for the a 
decay of natural tungsten isotope 180 W was observed for the first time. The measured half-life 
Ti/2 — l-lta4(stat)±0.3(syst)xl0 18 yr is close to the theoretical predictions, and it would be 
the most rare a decay ever observed in nature. More conservatively our result can be treated as 
the lower limit on half-life of 180 W: lim 7^ /2 ( 180 W) > 0.7xl0 18 yr at 90% C.L. Besides, new T 1/2 
bounds have been set for a decay of 182 W, 183 W, 184 W and 186 W at the level of (0.8- 1.8) xlO 20 
yr. All these limits are higher than those obtained in previous work [6] and are the most 
stringent bounds on Ty 2 for any a decaying nuclides. 

To confirm existence of a activity of 180 W, we are preparing measurements with other 
tungsten containing crystal scintillators: CdWC>4 (whose scintillation characteristics are better 
than those of currently used crystals), CaWC>4, and ZnWC>4. Observation of the a decay of 
180 W could be also checked with these crystals as bolometers [22, 38] and, apparently, sensitivity 
of such experiments would be enhanced by using crystals enriched/depleted in 180 W. 
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